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Abstract. The interpretation of microlensing results towards the Large Magellanic Cloud (LMC) still remains 
controversial. Whereas white dwarfs have been proposed to explain these results and, hence, to contribute signifi- 
cantly to the mass budget of our Galaxy, there are as well several constraints on the role played by white dwarfs. 
In this paper we analyze self-consistently and simultaneously four different results, namely, the local halo white 
dwarf luminosity function, the microlensing results reported by the MACHO team towards the LMC, the results 
of Hubble Deep Field (HDF) and the results of the EROS experiment, for several initial mass functions and halo 
ages. We find that the proposed log-normal initial mass functions do not contribute to solve the problem posed 
by the observed microlensing events and, moreover, they overproduce white dwarfs when compared to the results 
of the HDF and of the EROS survey. We also find that the contribution of hydrogen-rich white dwarfs to the 
dynamical mass of the halo of the Galaxy cannot be more than ~ 4%. 

Key words, stars: white dwarfs — stars: luminosity function, mass function — Galaxy: stellar content — Galaxy: 
dark matter — Galaxy: structure — Galaxy: halo 



1. Introduction 

White dwarfs are the most common remnants of stellar 
, evolution. Moreover, since white dwarfs are long-lived ob- 
' jects and the physical processes governing their evolution 
, are relatively well understood — at least up to moderately 
' low luminosities — they provide us with an invaluable tool 
, to study the evolution and structure of our Galaxy. In fact, 
the disk white dwarf luminosity function has become an 
important tool to determine some properties of the local 
neighborhood, such as its age (Winget et al. 1987; Garcfa- 
Berro et al. 1988; Hernanz et al. 1994), or the past history 
of the star formation rate (Noh & Scalo 1990; Di'az-Pinto 
et al. 1994; Isern et al. 1995a, b). This has been possi- 
ble because now we have improved observational lumi- 
nosity functions (Liebert, Dahn & Monet 1988; Oswalt et 
al. 1996; Leggett, Ruiz & Bergeron 1998) and because we 
have reliable cooling sequences — see, for instance, Salaris 
et al. (2000), and references therein. 

Although the situation for the disk white dwarf pop- 
ulation seems to be already clear and well understood, 
this is not the case for the halo white dwarf population. 
Moreover, the discovery of microlenses towards the Large 
Magellanic Cloud (Alcock et al. 1996, Alcock et al. 2000; 
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Lasserre et al. 2001) has generated a large controversy 
about the possibility that white dwarfs could be responsi- 
ble for these microlensing events and, thus, could provide a 
significant contribution to the mass budget of our Galactic 
halo. Nevertheless, white dwarfs as dark matter candi- 
dates are not free of problems, since an excess of this kind 
of objects would necessarily imply an overproduction of 
low-mass main sequence red dwarfs and high-mass stars 
that could eventually explode as Type II supernovae. In 
order to solve these problems Adams & Laughlin (1996) 
and Chabrier, Segretain & Mera (1996) proposed non- 
standard initial mass functions in which the formation of 
both low and high mass stars was suppressed. Besides the 
lack of evidences favoring such ad-hoc initial mass func- 
tions, they are not free of problems either. In particular, 
the formation of a typical {M ~ 0.6 Mq) white dwarf is 
accompanied by the injection into the interstellar medium 
of a sizeable amount of mass (~ 1.5 Mq on average). Since 
Type II supernovae are suppressed in ad-hoc initial mass 
functions, there is not enough energy to eject this mat- 
ter into the intergalactic medium and the ejected mass, 
which is significantly contaminated by metals (Abia et al. 
2001; Gibson & Mould 1997), cannot be accomodated into 
the Galaxy (Isern et al. 1998). Finally, an excess of white 
dwarfs also translates into an excess of binaries contain- 
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IMF 



Reference 



Standard 



CSMl 



CSM2 



AL 



' (M/Mo)-^-2 for M < 1.0 M© 
(M/Mo)-^ '^ for 1.0 Mo < M < 10.0 M® Scalo (1998) 
O.4(M/M0)-2-3 for M > 10.0 M© 



exp [-{M,ut/Mf^{M/M@f^] 
Mcut = 2.0 Mq, /3i = 2.2 y (32= 5.15 

exp [-{M,ut/Mf'{M/MQf^] 
Mcut = 2.7 Mq, /3i = 2.2 y /Ja = 5.75 



exp[l-^(log(f^J^] 
Mcut = 2.3 , y o- = 0.44 



Chabrier et al. (1996) 
Chabrier et al. (1996) 
Adams & Laughlin (1996) 



Table 1. The different IMFs used in this paper. 



ing such stars with the subsequent increase of Type la 
supernova rates which, ultimately, results in an increase 
in the abiuidances of the elements of the iron peak (Canal, 
Isern & Ruiz Lapuente 1997). All these arguments have 
forced the search for other possible explanations, such as 
self-lensing in the LMC (Wu 1994; Salati et al. 1999), or 
backgroiuid objects (Green & Jedamzik 2002) which have 
not been yet totally discarded. 

The suggestion of the MACHO team (Alcock et al. 
1997, 2000) that white dwarfs contribute significantly to 
the mass budget of the Galactic halo has motivated a 
large number of observational searches (Knox, Hawkins & 
Hambly 1999; Ibata et al. 1999; Oppcnhcimer et al. 2001; 
Majewski & Siegel 2002, Nelson et al. 2002) for these elu- 
sive white dwarfs. Also several theoretical works (Reyle, 
Robin & Crezc 2001; Flynn, Holopainen & Holmberg 
2003) have analyzed this possibility. However, the contro- 
versy of whether or not white dwarfs can provide a signifi- 
cant contribution to the Galactic dark matter is still open 
and deserves some more attention. In Isern et al. (1998) we 
analyzed the halo white dwarf population. In particular, 
we computed, assuming a standard initial mass function 
and updated models of white dwarf cooling, the expected 
luminosity function, both in luminosity and in visual mag- 
nitude, for different star formation rates. We showed that 
a deep enough survey (limiting magnitude <; 20) could 
provide important information about the halo age and the 
duration of the formation stage. We also showed that the 
number of white dwarfs produced using the proposed bi- 
ased IMFs could not represent a large fraction of the halo 
dark matter if they were constrained by the observed lu- 
minosity function of halo white dwarfs. However, within 
the approach adopted there the biases introduced by the 
sample selection procedures were not taken into account. 
More recently, we have analyzed (Torres et al. 2002) the 
sample of Oppenheimer et al. (2001). In this paper we ex- 
amine in detail the results of the MACHO team (Alcock 
et al. 2000) carefully taking into account the observational 



biases, thus updating our previous papers on this subject. 
We also study the number of white dwarfs which could be 
potentially found in the HDF (Ibata et al. 1999). Finally 
we also analyze the very recent observational results of the 
EROS team (Goldman et al. 2002), which set a very strin- 
gent upper limit to the white dwarf content of the Galactic 
halo. All these analyses are done by making use of a Monte 
Carlo simulator (Garcfa-Berro et al. 1999; Torres et al. 
1998). This is an important issue since white dwarf popu- 
lations are usually drawn from kinematically selected sam- 
ples (white dwarfs with relatively high proper motions). 
Therefore, some kinematical biases or distortions are ex- 
pected. A Monte Carlo simulation of a model population 
of white dwarfs is expected to allow the biases and effects 
of sample selection to be taken into account, so the prop- 
erties of the real sample could be corrected — or, at least, 
correctly interpreted — provided that a detailed simula- 
tion from the very early stages of source selection is per- 
formed accurately. Our paper is organized as follows. In 
section §2 we describe in full detail the Monte Carlo code. 
In §3 we present the results, whereas in §4 our conclusions 
are summarized. 

2. The model 

In this section we discuss the main ingredients of our 
Monte Carlo simulator. Since we want to self-consistently 
simulate simultaneously four different results, namely, the 
local halo white dwarf luminosity function, the microlens- 
ing results towards the LMC, the results of the HDF, and 
the results of the EROS team, and each one of these sim- 
ulations requires slightly different inputs, we will describe 
them in separate subsections. All of them, however, share 
some common ingredients like a random number genera- 
tor, which is always at the heart of any Monte Carlo sim- 
ulation. We have used a random number generator algo- 
rithm (James 1990) which provides a uniform probabihty 
density within the interval (0, 1) and ensures a repetition 
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period of ^ 10^*, which is enough for our purposes. Each 
one of the Monte Carlo simulations discussed in section 3 
below consists of an ensemble of 40 independent realiza- 
tions of the synthetic white dwarf population, for which 
the average of any observational quantity along with its 
corresponding standard deviation were computed. Here 
the standard deviation means the ensemble mean of the 
sample dispersions for a typical sample. 

2.1. The local halo white dwarf luminosity function 

We have considered a typical spherically symmetric halo. 
One of the most commonly used models of this type is 
the isothermal sphere. The density profile of the luminous 
halo is given by the law 



P{r) = Po- 



(1) 



where a ~ 5 kpc is the core radius, po is the local halo 
density and Rq =8.5 kpc is the galactocentric distance of 
the Sun. We randomly choose three numbers for the spher- 
ical coordinates (r, 6, </>) of each star of the sample within 
approximately 350 pc from the sun, according to Eq. (1). 
Afterwards we draw another pseudo-random number in 
order to obtain the main sequence mass of each star, ac- 
cording to one of the four model initial mass functions that 
will be studied here. Those initial mass functions are, re- 
spectively, a standard IMF (Scalo 1998), the biased IMF 
of Adams & Laughlin (1996) and the two ad-hoc IMFs 
proposed by Chabrier et al. (1996). All these mass func- 
tions are summarized in table 1. Once the mass of the pro- 
genitor of the white dwarf is known we randomly choose 
the time at which each star was born (ib). For this pur- 
pose we assume that the halo was formed 14 Gyr ago (see, 
however, section §3 below) in an intense burst of star for- 
mation of duration ^ 1 Gyr. Given the age of the halo, 

and the main sequence lifetime as a function of the 
mass in the main sequence (Iben & Laughlin 1989) we 
know which stars have had time enough to become white 
dwarfs, and given a set of cooling sequences (Salaris et al. 
2000) and the initial to final mass relationship (Iben & 
Laughlin 1989), which are their luminosities. 

The velocity distribution has been modeled according 
to a Gaussian law (Binney & Tremaine 1987): 



f{Vr,Vt) 



1 



1 



(27r)3/2 ara^ 



(2) 



where ct^ and at — the radial and the tangential veloc- 
ity dispersion, respectively — are related by the following 

expression: 



a. 



■ 



1 - 



2 , r dja^) 



(3) 



which, in a first approximation, leads to = at = 
— see, for instance, Binney & Tremaine (1987). For the 
calculations reported here we have adopted a circular ve- 
locity Vc = 220 km/s. Prom these velocities we obtain the 



heliocentric velocities by adding the velocity of the LSR 
^^LSR = —220 km/s and the peculiar velocity of the sun: 
(t/0, 1^0,1^0) = (10.0,15.0,8.0) km/s (Dehnen & Binney 
1998). Since white dwarfs usually do not have determina- 
tions of the radial component of the velocity, when needed 
for the observational comparison the radial velocity is 
eliminated. Moreover, we only consider stars with tangen- 
tial velocities in the range 250 ^vt^ 750 km/s. Stars with 
velocities smaller than 250 km/s would not be considered 
as halo members, whereas stars with velocities larger than 
750 km/s would have velocities exceeding 1.5 times the es- 
cape velocity, which we obtain from Binney & Tremaine 
(1987): 



v', =2V,'[1 + Hn/r)] 



(4) 



where r, ~ 41 kpc is the radius of the galactic halo. 

In order to build the white dwarf luminosity function 
using the l/Vma,^ method (Schmidt 1968) a smaller sample 
of white dwarfs must be culled from the original sample 
and in order to do this a set of restrictions in visual magni- 
tude and proper motion must be adopted. The restriction 
in magnitude will be discussed in §3 below. Regarding the 
proper motion cut we have chosen /i > 0.16" yr^^ as it 
was done in Oswalt et al. (1996) and in Garcia Berro et 
al. (1999). Besides, the l/Knax method requires that all 
the objects belonging to the restricted sample must have 
known parallaxes. This, in turn, means that all the white 
dwarfs belonging to this sample are within a sphere of 
radius of roughly 200 pc centered in the location of the 
sun. 



2.2. Microlensing events towards the LMC 

In order to produce a set of microlensing events towards 
the LMC we also need to simulate the characteristics 
of the white dwarf population towards the LMC. For 
this purpose we generate the three galactic coordinates 
{r,l,b) of the white dwarfs of the galactic halo inside a 
small pencil of 4° x 4° centered in the LMC location, 
{l,b) = (280°, -33°). The I and b distributions are prac- 
tically uniform in this small window. The radial coordi- 
nate is always smaller than the outer limit of the halo 
(r < 41 kpc) and according to the radial distribution 



P = Po{RQ/ry 



(5) 



where 7 = 3.4 — see, however, §3.7. We have chosen this 
distribution instead of that of equation (1) because in this 
way the number of microlensing events is maximized. In 
fact the mass distribution of microlenses does not neces- 
sarily follow the mass distribution of the luminous halo. 
However, in a first step (see §3.6 below), we have nor- 
malized the density of white dwarfs obtained from this 
distribution to the white dwarf density of the local neigh- 
borhood, n ~ 9.0 X 10"^ pc-3 for \og{L/LQ) Z -3.5 
(Tbrres et al. 1998). The velocity dispersions in this case 
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Fig. 1. Luminosity function of halo white dwarfs for several limiting magnitudes. 



are determined from Markovic & Sommer-Larsen (1996). 
For the radial velocity dispersion we have: 



— ^2,2 



11 r- 
arctan — - 

2 TT \ I 



(6) 



where ao = 80kms ^, a+ = 145kms ^, ro = 10.5 kpc 
and I = 5.5 kpc. The tangential dispersion is given by: 



1, 



+ 



r Aal 



where 
dr 



1 r cr| 
■^Tl + [(r-ro)/Z]2 



(7) 



(8) 



In order to decide which white dwarf could potentially 
produce a microlensing event a magnitude cut must be 
adopted. If the white dwarf is brighter than this limit- 
ing magnitude it could be potentially detected and, con- 
sequently, would not be a genuine microlensing event. 
Wc have explored a wide range of possibilities, namely 
m^"' = IT.S^'^s^ 22.5™'^s and 27.5'"''s, but, as a fiducial 
value, we have taken = 17. 5™*^, which is the value 
adopted by Alcock et al. (2000). We also need to simulate 
the population of stars of the LMC itself. In order to pro- 
duce such a population we distribute the number of mon- 
itored point sources (^ 1.2 x 10'') according to the LMC 
model of Gyuk et al. (2000). That is, we have adopted a 
scale lenght, i?d = 1-57 kpc, a scale height = 0.3 kpc. 



an inclination angle i = 30°, a distance L = 50 kpc, a 
position angle, (j) = 170°, and a tangential heliocentric ve- 
locity, vt = 336 kms~^. The next step is to determine if 
there exists a white dwarf in the line of sight to a given 
simulated star in the LMC. Wc consider a white dwarf to 
be responsible of a microlensing event if the angular dis- 
tance between the white dwarf and the monitored star is 
smaller than the Einstein radius 0e = Re/Dol, where Re 
is the Einstein radius and -Dql is the distance between the 
observer and the lens. This, of course, can happen at any 
time during the total monitoring period of 5.7 yr, due to 
the proper motion of the lenses, which is obtained from 
the previous equations. 

The optical depth is obtained, following Alcock et al. 
(2000), using the expression 



1 TT x ^ 



ti 



(9) 



where E = 6.12 x 10^ is the total exposure (in star- years), 
ti is the Einstein ring diameter crossing time, and e{ti) is 
the detection efficiency. The detection efficiency has been 
modelled as: 



e{i) = 



0.43 e- 



-(ln(t7T„))='-'58/0.87 



t>T„ 



0.43e-|in(*7T„,)|-='Vii.i6, i<T„ 



(10) 



where Tm = 250 days. This expression provides a good fit 
to the results of Alcock et al. (2000). 
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Fig. 2. Luminosity function of halo white dwarfs for several IMFs. 



2.3. The Hubble Deep Field simulation 



For the simulation of the HDF we have distributed stars 
in a window of 1.2° x 1.2° centered around = 
(125°, 55°). The radial distribution is, again, according 
to equation (5) within the outer halo limit. The veloci- 
ties are consequently drawn from equations (2) and (6) to 
(8). This simulation has been also normalized to the local 
density given by the halo simulation. 



2.4. The EROS simulation 



Finally, for the simulation of the EROS results we have 
distributed stars in a window of 190°^ in the Southern 
Galactic Hemisphere (-79° < 6gai < -48°), in the fol- 
lowing strips along the a coordinate, Ai5 = 1.4° wide: 
22iii6min ^ a < 3'^44™ at 6 = -44° 45', 23'^ 31™ < 
a < lh34min at 5 = -40° 09' and 22i>24™ < a < 
3^28™ at d = -38° 45', and 187°^ in the Northern 
Galactic Hemisphere (41° < 6gai < 59°) with lO'^ 57"^'" < 
a < 13^ 23""" at 6 = -12° and 10^ 57™"^ < a < 121^ 53™'" 
at ^ = —4° 36'. The radial distribution and the velocity 
distribution are, again, the same used for the HDF simula- 
tion. Also, the density of white dwarfs has been normalized 
to the local density given by the halo simulation. 



3. Results 

3.1. The halo white dwarf luminosity function 

One of the most serious problems that is found when de- 
termining the observational white dwarf luminosity func- 
tion of halo white dwarfs is that the real limiting mag- 
nitude used in these studies is highly uncertain. We have 
conducted a series of simulations with different limiting 
magnitudes to determine which would be the limiting 
magnitude able to reproduce the observational luminos- 
ity function. In figure 1 we show the white dwarf lumi- 
nosity functions obtained for several limiting magnitudes 
(16, 17, 18 and 19™''^) as solid lines. The luminosity func- 
tion of Torres et al. (1998) is shown as a dashed line for 
comparison purposes. Each panel is clearly marked with 
its corresponding limiting magnitude. The adopted halo 
age in all the cases was 14 Gyr. Also the adopted IMF is 
the standard one in all four simulations. The error bars of 
each luminosity bin were computed according to Liebert 
et al. (1988): the contribution of each star to the total er- 
ror budget in its bin is conservatively estimated to be the 
same amount that contributes to the resulting density; the 
partial contributions of each star in the bin are squared 
and then added, the final error being the square root of 
this value. This procedure is followed for each one of the 
40 realizations of the Monte Carlo simulation. After doing 
this the ensemble average of the dispersions is computed. 
Obviously the larger the magnitude limit the fainter the 
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Fig. 3. Microlensing optical depth towards the LMC as a function of the limiting magnitude for several IMFs. 



white dwarfs we detect. If we disregard as non significa- 
tive the bin in which we only detect on average one white 
dwarf we sec that the limiting magnitude that best re- 
produces the luminosity function of Torres et al. (1998) 
is my™ = 17™^s. Therefore for the rest of the simulations 
we adopt this value. It should be mentioned as well that 
in ordca' to detect the cut-off of the luminosity function a 
limiting magnitude of 25"^'*^ should be adopted. It is worth 
noting at this point that since the heliocentric velocity of 
halo white dwarfs is considerable, the proper motion cut 
plays a limited role. In fact the proper motion cut only 
affects the total number of white dwarfs in the sample 
but not the shape of the luminosity function. In contrast, 
as wc have seen, this is not the case for the cut in magni- 
tude. This is the same as to say that the proper motion cut 
equally affects all the luminosity bins. It is as well interest- 
ing to note here that the value of the limiting magnitude 
that we have found so far is in close agreement with the 
cut adopted by Alcock et al. (2000) for the monitoring of 
stars within the MACHO project, which is 17.5™'*^. 



3.2. The role of the IMF 

In figure 2 we show the luminosity functions obtained with 
a Salpeter-like mass function, the two log-normal IMFs 
proposed by Chabrier et al. (1996) — CSMl and CSM2, 
respectively — and the IMF of Adams & Laughlin (1996) 
— AL. As this figure clearly shows, the derived luminos- 



ity functions are not very sensitive to the precise shape of 
the IMF. Moreover, the completeness, as measured by the 
{V/Vynax) method, seems to be similar in all four cases. 
Only in the cases CSMl and, more apparently, in the 
CSM2 case there is a slight underproduction of luminous 
white dwarfs. However, the significance is only marginal 
and, therefore, in order to constrain the IMF of the galac- 
tic halo using intrinsically bright white dwarfs deeper sur- 
veys arc needed. We will come back to this issue when 
studying the HDF simulation in §3.4 below. 

As previously stated in §2.1, all the Imuinosity func- 
tions obtained here have been normalized to the local den- 
sity of halo white dwarfs obtained by Torres et al. (1998), 
n ~ 9.0 X 10^^ pc^'' for log(L/L0) > —3.5, which for a 
typical value of the mass of white dwarfs (~ 0.6 Mq) cor- 
responds to a density of halo baryonic matter in the form 
of white dwarfs of ~ 7.2 x 10~^ Mq pc~^. However, from 
our simulations we can derive which would be the total 
density of baryonic matter in the galactic halo within 300 
pc from the sun. We obtain po = 2.2 x 10~^ Mq pc~^ for 
the standard IMF. 2.0 x lO'^ Mq pc^^ for the CSMl case, 
1.1 X 10-2 Mq pc~'-^ for the CSM2, and 2.7x10-3 M© pc-^ 
for the AL simulation. These values, in turn, correspond 
to a fraction rj of baryonic dark matter of 0.03, 0.25, 1.40 
and 0.33, respectively. As it can be seen the differences 
between all the IMFs analyzed here are considerable. For 
instance, for the CSM2 case we would have more matter 
than needed, whereas the CSMl and AL mass functions 
would lock a sizeable fraction of dark matter in the form 
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Oil 
0.564 
0.016 
3.09 
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22.5 27.5 



Oil 
0.568 
0.014 
3.14 
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Oil 
0.556 
0.012 
5.05 
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75.8 89.3 



2i4 
0.590 
0.030 
1.53 
222 
36.3 



li3 
0.582 
0.024 
2.12 

243 
52.6 



1 i2 
0.601 
0.010 
6.19 

291 
70.8 



7i3 
0.598 
0.093 
0.51 

224 



6i4 
0.613 
0.035 
1.50 

251 



21.7 41.7 



2 i2 
0.595 
0.021 
3.04 
308 
70.9 



5i5 
0.622 
0.035 
1.63 

272 
28.8 



3i3 lil 

0.636 0.614 

0.030 0.012 

1.65 5.10 

232 293 



42.7 



69.2 



0.101 0.079 0.102 0.209 0.230 0.177 0.466 0.353 0.308 0.308 0.228 0.163 



Table 2. Summary of the results obtained for the simulation of microlenses towards the LMC for an age of the halo 
of 14 Gyr, different model IMFs, and several magnitude cuts. 



of main sequence stars, stellar remnants and in the cor- 
responding ejected mass. Finally the standard IMF only 
allows for a modest 3% of the required dark matter. In 
this case, moreover, roughly 1/3 of the stellar content cor- 
responds to white dwarfs. 

3.3. Microlenses towards the LMC 

In figure 3 we show the microlensing optical depth towards 
the LMC obtained in our simulations normalized to the 
value derived by Alcock et al. (2000), tq = 1.2 x 10""^, 
as a function of the magnitude cut for detection of white 
dwarfs (see section §3.1). As it can be seen in this figure 
the contribution of white dwarfs to the optical depth de- 
pends sensitively on the adopted IMF. For the standard 



IMF we derive a contribution of roughly 10%, no matter 
what the magnitude cut is, whereas for the biased IMFs 
we obtain contributions which are, typically, of 23% for 
the CSMl case, 35% for the CSM2 case, and 22% for the 
AL simulation. These values are relatively constant for 
large enough magnitude cuts and, hence, for realistic val- 
ues of the magnitude cut they can be considered as safe. 
Conversely, the microlensing optical depth is a robust in- 
dicator of the density of the microlenses. It is important 
to realize that none of the adopted IMFs is able to repro- 
duce satisfactorily the value found by the MACHO team, 
in spite of the very extreme assumptions adopted for de- 
riving the log-normal biased IMFs, which were especially 
taylored to reproduce the microlensing results. 
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Fig. 5. Distribution of the whole white dwarf population — solid lines — and of the white dwarfs responsible for 
microlensing events towards the LMC — dashed lines — as a function of their visual magnitude for several IMFs. 



Nevertheless the information that can be derived from 
our Monte Carlo simulations is far more complete. A sum- 
mary of the results can be found in table 2, where we show 
the number of microlensing events, the average mass of the 
microlenses, the average proper motion and distance, the 
average tangential velocnty of tlic^ microlenses, the corre- 
sponding Einstein crossing times and, finally, the contri- 
bution to the optical depth, all of them for three selected 
magnitude cuts. As it can be seen in the table 2 none 
of the IMFs is able to reproduce the observed number of 
microlensing events (13 to 17, depending on the selection 
criteria) found by Alcock et al. (2000). Even in the case 
of the CSM2 simulation, which corresponds to the most 
extreme assumption on the IMF, we obtain only 7 ± 3 mi- 
crolensing events in the best of the cases. This explains 
why the contribution to the microlensing optical depth of 
the white dwarf population is within the range 10% to 
35%, depending on the adopted IMF, suggesting that a 
sizeable fraction of the microlensing events could be ei- 
ther self-lensing in the LMC (Salati et al. 1999) or due to 
background objects (Green & Jedamzik 2002). 

Another important information which can be readily 
obtained from the simulations presented so far is the dis- 
tribution of Einstein crossing times. Such distributions for 
the 40 independent realizations are shown in figure 4 for 
the four cases studied here as solid lines. We have chosen 
the brighter of our magnitude cuts in order to allow for 
a larger number of events. Note, however, that this mag- 



nitude cut is exactly the same adopted by the MACHO 
team. Also shown, as dashed lines, are the distribution of 
Einstein crossing times obtained by the MACHO team. All 
the distributions have been normalized to unit area and, 
hence, are frequency distributions. The average Einstein 
crossing times for each of the simulations can be also foimd 
in table 2. As it can be seen in figure 4 the MACHO 
team detects microlensing events with larger durations 
than those obtained in the simulations of the log-normal 
IMFs of Chabrier et al. (1996) and Adams & Laughlin 
(1996). Although the statistics is poor, we can compare 
the average crossing times obtained here with the average 
crossing time obtained by the MACHO team which is of 
76 ±23 days. These average crossing times are respectively 
^ 61, 36, 22 and 29 days. In all the cases the typical stan- 
dard deviation is of about 15 days. We emphasize here 
that these averages are the result of an ensemble of forty 
simulations. Clearly, biased mass functions yield consider- 
ably smaller average Einstein crossing times than those of 
the MACHO team, whereas a standard IMF yields a value 
which is actually within the la error bars of the observed 
one. Since (tE)^ is an indicator of the mass of the lens it 
follows that the average mass of the lenses for biased mass 
functions must be different of that of the standard one. 
And this is, indeed, the case as it can be seen in table 2. 
This behavior is not surprising since the three log-normal 
mass functions partially inhibit the formation of low mass 
white dwarfs. In all these three cases the observed average 
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AL 




Age 


12 Gyr 


14 Gyr 


16 Gyr 


12 Gyr 


14 Gyr 


16 Gyr 


(iVwD) 


8±3 


7±3 


7±3 


5±5 


5±5 


3±3 


(m) (M/Mo) 

(m> ("yr-^) 


0.620 


0.598 


0.589 


0.642 


0.622 


0.650 


0.039 


0.093 


0.041 


0.043 


0.035 


0.026 


(d) (kpc) 


1.24 


0.51 


1.04 


1.07 


1.63 


1.66 


(Vtan) (kms-1) 


230 


224 


202 


219 


272 


203 


(iE) (d) 


27.0 


21.7 


20.9 


38.1 


28.8 


48.4 


(r/^o) 


0.491 


0.466 


0.415 


0.323 


0.308 


0.310 



Table 3. Summary results of the microlensing events towards the LMC for a limiting magnitude of my = 17.5' 
and different halo ages and IMFs. 



Einstein crossing time is beyond the la error bars 15 
days). Hence, although this mass functions yield a signif- 
icantly higher number of microlensing events they have 
considerably smaller Einstein crossing times. 

Another interesting fact which results from a careful 
study of table 2 has to do with the tangential velocities of 
the simulated microlensing events. As it can be seen the 
average tangential velocity of the microlenses is typical of 
the halo population, with values close to the canonical one 
of 220 km/s. But, on the contrary, the average distances 
of the microlenses do depend on the adopted IMF. For 
instance, the microlenses produced by the standard IMF 
are located, on average, at larger distances (3.14 kpc), 
than those of the biased IMFs (2.12, 1.50 and 1.65 kpc, 
respectively). In fact, the larger the adopted mass cut of 
the IMF is (see table 1), the closer the microlenses are. 
This, in turn, translates directly into the average proper 
motions of the lenses, since, as already mentioned, the 
tangential velocities of the lenses do not differ much in 
the simulations reported here. Finally, as expected, the 
average mass of the microlenses differs considerably in all 
the cases studied here. For the standard IMF we obtain 
a mass of 0.56 Mq, close to the canonical average mass of 
disk white dwarfs. For the biased IMFs we obtain values 
of 0.58 Mq, 0.61 Mq and 0.64 Mq, respectively, which re- 
flect, as already pointed out, the adopted mass cut for the 
IMF. 

In figure 5 we show the relative distribution (normal- 
ized to unit area) in magnitudes of the whole population of 
white dwarfs — solid lines — and of those white dwarfs re- 
sponsible of the microlensing events — dashed lines. Also 
shown in each panel is the total number of white dwarfs 
found in a typical Monte Carlo realization. Again, the dis- 
tributions have been normalized to unit area. As it can be 
seen in this figure, most of the white dwarfs have apparent 



magnitudes brighter than 25™**^ and, hence, the magni- 
tude cut only plays a significant role in the case in which 
the magnitude cut is 21 .'t>™^^^ which is not very realistic. 
Therefore the results quoted in table 2, do not depend 
much on the adopted magnitude cut. 

Finally, we have analyzed the dependence of these re- 
sults on the adopted age of the halo. In order to do that 
we have chosen halo ages of 12, 14, and 16 Gyr. In all 
these simulations we have adopted the same magnitude 
cut my"' = 17.5™*^s, xhe results obtained in these sets of 
simulations are shown in table 3. From an analysis of this 
table we see that the results do not depend much on the 
adopted halo age. In particular the number of microlensing 
events remains almost constant regardless of the precise 
value of the age of the halo. Also, and most importantly, 
the same behavior can be observed for the computed mi- 
crolensing optical depth. Thus, the actual values of the 
number of microlensing events and of the microlensing op- 
tical depth derived here are quite robust. 

3.4. The Hubble Deep Field 

The deepest optical images obtained up to now are those 
of the Hubble Deep Field. In spite of the very small area 
surveyed by the HDF (~ 4.4 arc min^), the limiting mag- 
nitude of my™ ~ 28™'*^ allows to probe a large volume. 
Therefore, these results complement those of the MACHO 
team and can provide us with relevant (and very valuable) 
information about the halo of the Galaxy. Nevertheless, 
the results obtained so far by different authors are not 
concluding. For instance Flynn et al. (1996) studied the 
HDF north and were able to distinguish between stars and 
galaxies down to magnitudes as faint as / = 26. Their 
selection criteria in the search of white dwarfs are sum- 
marized in the first row of table 4. They did not find any 
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Fig. 6. Theoretical color-magnitude diagram of a typical Monte Carlo realization of the HDF for several IMFs. 



object with V — I> 1.8 although their results were consis- 
tent with an upper limit of 3 white dwarfs. Later Mendez 
et al. (1996) detected 6 faint objects {V > 25) with color 
indexes in the range —0.5 < B — V < 0.5 or, equiva- 
lently, < V — I < 1.2. However, these objects could be 
as well non resolved galaxies. Finally, the most plausible 
analysis of the HDF is that of Ibata ct al. (1999). These 
authors determined the proper motions of the objects in 
the HDF from two epoch observations. The time baseline 
was 2 yr. In this way they were able to discriminate be- 
tween nearby objects and galaxies. The maximum number 
of white dwarfs found by Ibata et al. (1999) is 4. However, 
one of these identifications corresponds to a re-discovery 
of a previously known white dwarf, whereas the rest of the 
candidates still need spectroscopic follow-up observations. 
It is, thus, interesting to study which would be the pre- 
dictions of our Monte Carlo simulator for different IMFs. 

In figure 6 the theoretical color-magnitude diagram of 
a typical Monte Carlo realization of the HDF is shown. As 
it can be seen in this figure most white dwarfs populate the 
coolest and reddest portion of the cooling isochrone. This 
is a consequence of the fact that the characteristic cooling 
times increase considerably at low luminosities due to both 
the release of latent heat upon crystallization and of grav- 
itational energy due to carbon-oxygen separation (Isern 
et al. 2000). The characteristic "z" -shaped feature at the 
lowest luminosities is due to the contribution of massive 
white dwarfs and not to the blue turn reported by Hansen 
(1999) — see Isern et al. (1998) and Salaris et al. (2000) 



Author 



Selection Criteria 



N, 



Obj 



Flynn et al. (1996) 



Mendez et al. (1996) 



Ibata et al. (1999) 



24.63 < 7 < 26.30 
V-I>1.8 
25 < < 30 
-0.5<B-F<0.5 
Q<V-I <1.2 
27 < F < 28 
-0.2 <V -I <1.0 



< 3 



Table 4. Summary of the HDF results. The authors, the 
selection criteria and number of objects are indicated. 



for a comprehensive explanation of this feature — which 
our cooling sequences reproduce well but for larger ages 
of the Galactic halo. Note as well that for biased mass 
functions the contribution of low mass stars has been sup- 
pressed. This, in turn, implies that these mass functions 
should provide a very small number of white dwarfs with 
moderately high luminosities (Isern et al. 1998) and, in- 
deed, this is the case. Of course for larger mass cuts the 
number of white dwarfs on the hot portion of the cooling 
isochrone should be smaller and, hence, this is the reason 
why in the AL case and the CSM2 simidation white dwarfs 
tend to concentrate at the hook of the cooling isochrone. 

In figure 7 a typical realization of the simulation of 
the HDF is shown. In this figure the box represents the 
observational selection criteria of Ibata et al. (1999), as 
shown in table 4. Since for biased mass functions most 
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Fig. 7. Simulation of the HDF. The box represents the observational selection criteria of Ibata et al. (1999). 



white dwarfs tend to be concentrated at the hook of the 
cooling isochrone and they are located at very different 
distances the result is that for these cases white dwarfs 
are distributed along an almost vertical strip, which is lo- 
cated at the position of the hook, V — I ^ 0.6. This could 
be a potential way to test the IMF of the Galactic halo 
in addition to the total number of white dwarfs found in 
the field, since, as expected, biased mass functions pro- 
vide a considerably larger number of white dwarfs when 
compared with the standard case. 

The fraction of dark matter of a given population is 
usually computed according to: 



J. _ 3A^objM j^q-0.6(tov-Mv+5) 



(11) 



where A^obj is the number of objects in the sample, pdm = 
8.0 X 10~'^pc~'^ is the local density of halo dark matter, 
M is the average mass of the white dwarfs in the HDF for 
a given simulation, my is the apparent magnitude cut of 
Ibata et al. (1999) as given in table 4, My is the absolute 
magnitude of the dimmest white dwarf in the simulated 
sample, = 4.4 arc min^ is the area surveyed by the 
HDF, and e is the fraction of white dwarfs within the 
selection criteria over the total number of white dwarfs. 
This fraction has been estimated to be e ~ 0.42 (Ibata 
et al. 1999), but its real value is uncertain. Consequently, 
this expression does not fully take into account all the 
biases introduced by the selection process and, hence, the 
adopte value of e can only be regarded as a relatively bona 



fide estimate. Instead, the dark matter fraction has been 
computed in the following way. 

The expected number of objects in the sample, is given 
by -^'obj = nVcS, where VcS is the effective volume of the 
sample and n = fpo/M. As shown in Eq. (11) the usual 
way to compute the effective volume of the sample is Ves = 
e^il.(P where d is taken to be logcZ = 0.2(mv — My + 5). 
Another, more accurate, way of computing the effective 
volume is to use the 1/Knax method as it follows. For each 
star of the sample we determine the maximum distance 
over which any star can contribute to the sample. 



W = min [7r-i(/z/^i); tt-HoO-^^™'"™) 
and the minimum distance as well: 



= max 



TT-i (/.//.„); ;r-il0°-2(™-™) 



where it is the; stcillar parallax, /i its proper motion, m 
the apparent magnitude, /xi and fiu are the low and high 
proper motion limits, if any, and mi and are the corre- 
sponding magnitude limit. The maximum volume in which 
a star can contribute is then 



3 

— (r 

^ V max 



and the number density of white dwarfs will be 
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IMF 


A^obj 


Kff (pc-^) 




e 


d (kpc) 


/ 


Standard 


2±2 


6360 


0.396 


0.69 


3.7 


0.04 


CSMl 


4±4 


2051 


0.466 


0.79 


2.5 


0.20 


CSM2 


30 ±6 


874 


0.475 


0.68 


1.9 


4.16 


AL 


3±3 


2180 


0.531 


0.80 


2.6 


0.14 



Table 5. Summary of the results for the HDF obtained with our Monte Carlo simulator for the different proposed 
IMFs. 



The effective volume of the sample can be then computed 
as: 



obj 



Veff = 

n 

The value of {V/Vmax) is a measure of the completeness 
of the sample, which is related to e by 



V 



0.5 



We remind here that for a complete and homogeneous 
sample {V/Vmax) =0.5. 

The results for the HDF simulations are also shown 
in table 5. For each of the cases we show the number of 
expected white dwarfs in the HDF, A'obj, with its cor- 
responding standard deviation, the effective volume sur- 
veyed, the value of {V/Vmax), the estimated completeness, 
the average distance of the white dwarfs in the field, and 
the corresponding fraction of dark matter derived from 
the simulations. An analysis of table 5 reveals that the 
expected number of objects is roughly 3 in all the cases 
except for the case CSM2, for which the value of expected 
white dwarfs is much larger. It is interesting to compare 
these values with the expected number of objects if all 
the dark matter of the Galaxy were in the form of white 
dwarfs. This number ranges, depending on the selection 
criteria in colors and magnitudes, from 9 to 12. Clearly 
even biased mass functions such as the AL and CSMl 
mass functions are not able to fill the halo with white 
dwarfs. Moreover, for the CSM2 case even if we adopt the 
lower limit of 24 objects in the field the result is clearly 
much larger than that needed to fill all the halo with white 
dwarfs. 

On the other hand, the effective volume surveyed by 
the HDF turns out to be dependent of the adopted model 
IMF. This would not be the case if Eq. (11) would have 
been adopted since in this case the effective volume sur- 
veyed only depends on the absolute magnitude of the 
faintest white dwarf in all the simulated samples which 
turns out to be ~ 16. 8™^^ — in accordance with the value 
adopted by Richer et al. (2000) and on the adopted 
magnitude cut (my = 28. 5™*^). It follows then that if 
this was the case, the radius of the effective volume would 



be in all cases ~ 2.2 kpc. The average proper motions 
obtained are in all four cases very difficult to measure: 
~ 26 mas yr~^, for the standard IMF, ~ 30 mas yr~-^ 
for the CSMl case, ~ 34 mas yr~i for the CSM2 mass 
function and ~ 32 mas yr~^ for the AL case. The typi- 
cal standard deviation is ±15 mas yr~^. It is important 
to realize that these values are in agreement with the ob- 
served value of 24 ± 10 mas yr'^ (Ibata et al. 1999). The 
completeness of the sample is roughly the same in all four 
cases, whereas for the average distance of the microlenses 
we have significant variations. In particular for the stan- 
dard IMF the average distance is considerably larger than 
in the rest of the simulations, being the case CSM2 an ex- 
treme case. Finally, the fraction of baryonic dark matter 
in the form of white dwarfs for the CSM2 case is totally 
incompatible with the observations, whereas in the rest of 
the cases is significantly smaller. For instance, for the case 
the CSMl and AL mass fimctions the fractions are com- 
parable, and for the standard IMF the fraction of baryonic 
dark matter only amounts to a modest 4%. Finally it is 
worth noticing here that if we had used the selection cri- 
teria of Flynn et al. (1996) we would have not found any 
white dwarf in the HDF, whereas if we have used those 
of Mendez et al. (1996) we would have found 1 ± 1 white 
dwarfs in HDF, almost independently of the adopted IMF. 
Hence, it seems that a significant fraction of the objects 
found by Mendez et al. (1996) should not be white dwarfs. 



3.5. The EROS results 

The EROS team has reported very recently (Goldman 
et al. 2002) the results of a proper motion survey which 
aimed to discover faint halo white dwarfs with high proper 
motions. However, they did not detect any candidate halo 
white dwarf even if the survey was sensitive down to 
y w 21 and to proper motions as large as ^ = 6.0" yr~^. 
Moreover, they found that the halo white dwarf contri- 
bution cannot exceed 5% at the 95% confidence level for 
objects with color index 1.0 < V — I < 1.5. It should 
be noted however that this last result is dependent on 
the adopted model of the Galaxy. In particular Goldman 
et al. (2002) adopted the biased mass function CSMl 
and in order to compare with the observational results 
they only simulated white dwarfs with magnitudes in 
the range 16.5 < My < 18.0 and color index within 
—0.5 < V — I < 1.5. It is therefore necessary to extend the 
study to other mass functions and, moreover, to the full 
range of white dwarf magnitudes and colors. Additionally, 



E. Garci'a-Berro, S. Torres & J. Isern: The halo white dwarf population 



13 




15 20 25 30 15 20 25 30 
V V 



Fig. 8. Simulation of the EROS results. See text for a detailed description. 



IMF 




Kff (pc^) 


{V/K.ax> 


e 


d (pc) 


/ 


Standard 


2± 1 


46577 


0.210 


0.42 


122 


0.01 


CSMl 


6±2 


13256 


0.115 


0.32 


81 


0.16 


CSM2 


12 ±4 


6042 


0.150 


0.30 


62 


0.55 


AL 


6±4 


10238 


0.140 


0.28 


74 


0.17 



Table 6. Summary of the results for the EROS experiment obtained with our Monte Carlo simulator for the different 
proposed IMFs. 



since the results of the EROS team are closely connected 

with those of the HDF studied in §3.4 before, it is inter- 
esting to study the results obtained with our Monte Carlo 
simulator. 

In figure 8 a typical Monte Carlo realization of the 
EROS results is shown, whereas in table 6 the average 
values for the ensemble of 40 independent realizations are 
also shown. The entries in table 6 are the same of table 
5, except for the distance which is expressed in pc. The 
selection criteria of the EROS collaboration are shown as 
dashed lines in figure 8. First, halo white dwarfs should 
have a reduced proper motion Hy = V+5 log7r+5 > 22.5, 
this restriction is shown in figure 8 as a dashed diag- 
onal line. Additionally they required 6.0" yr~^ > n > 
0.7" yr~^. Both limits are shown as well in figure 8 as 
a horizontal dashed lines. As it can be seen in table 6 
and in figure 8, in the region where the EROS experiment 
conducted their search for halo white dwarfs very few of 
them would be eventually found. In particular for the case 



in which a standard mass function is used only 2 of them 

(in the best of the cases) would be foimd. This is not 
the case for the biased mass functions CSMl and AL for 
which 6 white dwarfs could be presumably found, whereas 
for the CSM2 case this number increases to 12. Since the 
EROS experiment found none the only mass function that 
seems to fit the observational data is the standard one. In 
summary, both the results of the HDF, and of the EROS 
experiment point towards the same conclusion: that the 
IMF adopted for the halo should be the standard one. 

3.6. What if...? 

As it can be seen from tables 2 and 3 our results are com- 
patible with a halo where roughly 10% of the dark matter 
is in the form of white dwarfs. It is thus interesting to ask 
ourselves what would eventually be the observed luminos- 
ity function of halo white dwarfs if all the microlensing 
events (and the associated optical depth) are due to halo 
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Magnitude 

(iVwD) 

(m) (M/Mo) 

(m) ("yr-^) 

{d) (kpc) 

(Ftan) (kms-^) 

(te) (d) 



-^^gmag 22.5™^^ 27.5™''* 



Oil 
0.685 
0.158 
0.33 
248 
18.5 
0.064 



Oil 
0.660 

0.062 
1.21 
360 
23.2 

0.061 



Oil 
0.651 
0.019 
4.25 
386 
40.5 
0.060 



Table 7. Summary of the results obtained for the mi- 
crolensing events towards the LMC when the density pro- 
file of Navarro et al. (1997) is adopted. See text for more 
details. 



in the LMC and SMC are binary and, hence, it follows 
that a substantial fraction of single events must also be in 
the LMC or the SMC as first pointed out by Sahu (1994). 
However, it is true as well that this simulation reinforces 
the result that most of the microlensing events reported 
by the MACHO team arc not halo white dwarfs and that 
another explanation must be found for these results. 



white dwarfs. To this regard we have conducted a series 
of simulations in which we have produced as many white 
dwarfs as needed in order to reproduce the optical depth 
found by the MACHO team, which roughly corresponds 
to 13±1 microlenses. In all these simulations the standard 
IMF was adopted, since the biased mass functions AL and 
CSMl studied in the previous sections produce roughly 
the same results with regard to the white dwarf luminosity 
function. It is worth noticing as well that the ad-hoc mass 
functions studied here overproduce white dwarfs either in 
the HDF or in the EROS survey in sharp contrast with 
the observations. The white dwarf luminosity function was 
then computed using the restrictions for the visual mag- 
nitude and proper motion detailed in §2.1. As was done 
in §3.1, forty independent realizations of this sample were 
computed in order to obtain the average white dwarf lu- 
minosity function along with its corresponding standard 
deviation for each luminosity bin. The result is shown in 
figure 9 as a solid line. For the sake of the comparison the 
observational disk and halo white dwarf luminosity func- 
tions of Bergeron, Legget & Ruiz (1998) for the disk and 
of Torres et al. (1998) for the halo are also shown as dotted 
lines. As it can be seen in this figure, should white dwarfs 
be responsible of all the microlensing events the halo white 
dwarf luminosity function would be much closer to that 
of the disk rather than to that of the halo. It is, how- 
ever, worth mentioning at this point that this is indeed 
an overestimate since for at least 4 lensing events it has 
been already shown that the lenses are not in the Galactic 
halo. In particular, one of the lenses resides in the galac- 
tic disk (Alcock et al. 2001a). Moreover, three of these 
are binary events (Afonso et al. 1988; Bennet et al. 1996; 
Alcock et al. 2001b) in the LMC or SMC. It is obvious, 
that it is not possible that all lensing events with lenses 



3. 7. The role of the density profile 

In order to check the sensitivity of our results with respect 
to a different density profile we have done a final series of 
simulations in which the density profile of Navarro, Frenck 
& White (1997) was adopted. This density profile is char- 
acterized by the following expression 



p(r) 



Pcrit'^c 



{r/n){l + r/r,f 



(12) 



where Tg is the scale radius, 5c is an adimensional charac- 
teristic density, and pcrit = 3H^/8ttG is the critical den- 
sity of an Einstein-De Sitter universe. On the other hand, 
Ts is usually defined through the concentration parame- 
ter c = Ram/fs as a function of the virial radius, i?200- 
In order to correctly reproduce the rotation curve of our 
Galaxy the following set of parameters is usually adopted: 
c = 10, i?200 = 180 kpc, which results in fg = 18 kpc. In 
this simulation, and as discussed previously, the adopted 
initial mass function was the standard one. We only com- 
puted the number of microlenses towards the LMC and 
not the results of the HDF or of the EROS experiment. 

The different density profiles used here are displayed 
in figure 10. The solid line corresponds to the classical 
isothermal sphere density profile, the short dashed line to 
a power law with exponent 7 = 2.0, the long dashed line 
to a power law with 7 = 3.4 and the dashed-dotted line 
to the Navarro et al. (1997) density profile. The top panel 
shows the density profile as a function of the galactocentric 
r coordinate, whereas the bottom panel shows the density 
profile projected towards the direction of the LMC. As it 
can be seen in the top panel of figure 10, all the density 
profiles, including that of Navarro et al. (1997), are much 
more concentrated towards the center of the Galaxy than 
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Fig. 10. The density profiles used in the simulations. The 

top panel shows the density profiles from the center of 
the Galaxy, whereas the bottom panel shows the density 
profiles towards the direction of the LMC. 



the isothermal sphere. Additionally, the power law with 

7 = 3.4 decreases much faster than the other three den- 
sity profiles considered here. However, when projected to- 
wards the LMC, and for the distances of interest (those 
relevant for the study of the microlenses towards the LMC, 
obtained in the previous sections, 3 kpc), all four density 
profiles yield similar densities. Consequently, we do not 
expect large departures from what has been found in pre- 
vious sections, and this is indeed the case. In particular, 
the results obtained using the Navarro et al. (1997) den- 
sity profile arc shown in table 7, and it is evident from a 
cursory inspection that the results do not differ much from 
those shown in table 3. Hence, the same conclusions ap- 
ply, namely that the contribution of halo white dwarfs to 
the microlensing events towards the LMC must be small, 
independently of the adopted density profile. 

3.8. Miscellanea 

In order to check the robustness of our results we have con- 
ducted a series of numerical experiments in which, firstly, 
the distance to the LMC was varied by a 10% and, sec- 
ondly, the tangential velocity of the resulting white dwarf 
population was also varied by a 10%. In all the cases a 
standard IMF was used. The results are shown in table 8. 
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0.063 


0.101 


0.086 


0.092 


0.132 



Table 8. Summary of the results obtained for a limiting 

magnitude of rny = 17.5™^^ and difcrent values of the 
LMC distance and the tangential velocity 



Our fiducial model is shown in column 3. The second col- 
umn corresponds to a model in which a distance of 45 kpc 
was adopted, whereas the model in which a distance of 55 
kpc was used is shown in column 4. Column 5 and 6 corre- 
spond, respectively, to the models in which the tangential 
velocities were decreased (V^an) increased (Fj^n) ar- 
tificially by a 10%. As can be seen there, the differences 
in all the cases are relatively small. For instance decreas- 
ing the distance to the LMC from 50 to 45 decreases the 
contribution of white dwarfs to the observed microlensing 
optical depth by 6%, whereas the model in which the con- 
tribution of white dwarfs to the microlensing optical depth 
would be largest is that in which the tangential velocities 
of all the white dwarfs have been increased by 10%. In 
this case the contribution amounts to a 13%. Given all 
these considerations it seems unlikely that we could have 
underestimated the contribution of white dwarfs to the 
observed microlensing optical depth. 

4. Discussion and conclusions 

Using a Monte Carlo code we have computed self- 
consistently and simultaneously the theoretical expecta- 
tions of the number of white dwarfs for four different 
observational results, namely the reported microlensing 
events towards the Large Magellanic Cloud (Alcock et al. 
1997, 2000), the local halo white dwarf luminosity function 
(Torres et al. 1998), the results of the Hubble Deep Field 
(Ibata et al. 1999) and the results of the EROS experiment 
(Goldman et al. 2002). Our Monte Carlo simulator takes 
into account all the known observational biases and uses 
a thorough description of the properties of the Galactic 
halo. In our calculations we have simulated the halo white 
dwarf population for several halo ages (12, 14 and 16 Gyr). 
We have also used four different initial mass functions: a 
standard IMF (Scalo 1998), the two ad-hoc initial mass 
functions of Chabrier et al. (1996), and the biased initial 
mass function of Adams & Laughlin (1996) — see table 1. 

We have found that none but the standard IMF is able 
to reproduce simultaneously the observations of the HDF, 
of the EROS experiment and of the local halo white dwarf 
luminosity function. In particular, the most extreme initial 
mass function of Chabrier et al. (1996) largely overpro- 
duces white dwarfs for both the EROS survey and for the 
HDF. More precisely, the EROS experiment in their sur- 
vey did not find any candidate white dwarf whereas our 
simulations show that if this was the initial mass func- 
tion of the halo they should have found around 12 white 
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dwarfs. The same holds for the HDF, where a maximum 
of 4 white dwarfs should have been found, whereas our 
simulations using this IMF yield about 30 white dwarfs. 
Thus, this IMF is completely ruled out by observations. 
Regarding the more conservative IMF of Chabrier et al. 
(1996) our simulations predict 4± 4 white dwarfs for the 
HDF and 6 ± 2 for the EROS survey. Therefore, although 
this IMF is marginally compatible with the results of the 
HDF and with the local white dwarf luminosity function, 
it turns out that is incompatible with the results of the 
EROS survey. Since the Adams & Laughlin (1996) initial 
mass function is very similar to the Chabrier et al. (1996) 
mass function the same conclusions apply to this IMF. 

Regarding the microlensing events towards the LMC 
we have found that white dwarfs at most could be respon- 
sible of ~ 1 event for the case in which the standard IMF 
is adopted. This number increases up to 5 events for 
the cases in which the "moderate" log-normal mass func- 
tion of Chabrier et al. (1996) and of Adams & Laughlin 
(1996) are used. The contribution of these events to the 
observed optical depth tq is, respectively, t/to ^ 0.1, 0.2 
and 0.3. Hence, even if ad-hoc initial mass functions are 
adopted the total number of microlensing events produced 
by halo white dwarfs and the corresponding optical depth 
found by the MACHO team cannot be recovered when the 
observational biases and other constraints are taken into 
account. Moreover, the average duration of the simulated 
microlensing events for all the samples in which biased 
mass functions were used is considerably smaller than the 
observed one. In summary, the ad-hoc initial mass func- 
tions of Chabrier et al. (1996) and of Adams & Laughlin 
(1996), which were tailored to fit the observations of the 
MACHO team, fail to reproduce the rest of the observa- 
tional constraints and do not yield the required number of 
microlensing events (and, thus, do not significantly con- 
tribute to the optical depth). Consequently they should 
be disregarded until new (and unlikely) evidence could 
eventually appear. 

Wo would like to emphasize at this point that all the 
above mentioned results are almost completely indepen- 
dent of the adopted halo age. Moreover, in order to check 
whether or not our results are sensitive to the density pro- 
file of the halo, we have extended our simulations using the 
density profile of Navarro et al. (1997). We have found that 
the microlensing results are not sensitive to the choice of 
the density profile because at the distances of interest all 
the density profiles studied here are coincident. 

In summary, our results strongly suggest that a size- 
able fraction of the dynamical mass in the halo of our 
Galaxy cannot be in the form of old liydrogon-rich white 
dwarfs. More specifically, for the case in which a standard 
initial mass function is adopted this fraction turns out to 
be of the order of 4% if we trust the results of the HDF or 
as low as 1% if the results of the EROS team are adopted. 
In any case the percentage obtained here is well below 
the claim of the MACHO team that 20% of the dark halo 
of the Galaxy is tied up in half-solar mass stellar bodies 
(Alcock et al. 2000), since the constraints set by red star 



counts (Flynn et al. 1996) discard the possibility of low- 
mass main sequence stars. On the other hand our results 
arc in excellent agreement not only with the observations 
of the HDF and of the EROS experiment but also with 
the theoretical estimates of Graff et al. (1998), ^ 10%, 
Flynn et al. (2001), ^ 2%, and more recently Goldman et 
al. (2002), ^ 4%, and Zhao (2002), ^ 4%. Of course, the 
real nature of the reported microlensing events remains 
still a mistery and, hence, the last word on this issue has 
not been told yet. 
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